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Abstract 

The universal connected correlations proposed recently between eigenvalues 
of unitary random matrices is examined numerically. We perform an ensemble 
average by the Monte Carlo sampling. Although density of eigenvalues and a 
bare correlation of the eigenvalues are not universal, the connected correlation 
shows a universal behavior after smoothing. 
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Brezin and Zee |T]-§[ have discovered recently the universal behavior of the eigenvalue 
correlation of random matrices. Let <p be N x N hermitian matrices and let ith eigenvalue 
of cf) be denoted by Xi, pi or z/j. We consider the probability distribution with a weight 

P((f)) oc exp(-NTrV((j))) (1) 

where V{<fi) is an even polynomial of 0. The density of the eigenvalues p(X) and the corre- 
lation of the eigenvalues p(p, v) are defined as 

N 



p(A) = d ( 2 ) 

i=i 



^ N ^ N 

Pit*,") = (jr 53^-^)^53^-^)) (3) 

1=1 j=l 

where (. . .) is an ensemble average. Quantities p(X) and p(p, v) are not universal and depend 
on the details of P{4>), i-e. V((f)). Using orthogonal polynomials, one can calculate 
the correlation of the eigenvalues of random matrices. Using an ansatz for a form of the 
orthogonal polynomials, Brezin and Zee [U have discussed a connected part of p(p, u), which 
is defined by 

Pc(p, v) = p(p, v) - p{p)p{v) (4) 

This p c {p, v) oscillates wildly as a function of p and v. They have shown that, in the large 
N limit, the smoothed connected correlation p^ mooth (/z, v) has an universal form 

P s r oth (p, v) = -^v /(///a> u/a) > (5) 

with 



where ±a is the endpoint of the spectra of the eigenvalues. The smoothed connected cor- 
relation p™ ooth (yU, v) must be averaged over intervals Sp and 5u much less than 0(1) but 
larger than 0(1/N). Thus (|5]) is valid where \p— v\,\±a — p\ and | ± a — u\ are larger than 
0(1/N). The smoothed connected correlation p s c mooth (p, u) is universal in the sense that the 



function f(x,y) does not depend on the potential V at all. This universality has since been 
verified and derived by Beenakker H and by Eynard || using other methods. 

In this paper we study numerically the correlation of the eigenvalues to check the validity 
of the universality in equations (|^) and In our calculations, we limit the form of V(4>) 
to 

V{cj>) = v 2 <p 2 + + (7) 

We choose the matrix size N to be 100 and performed the ensemble average by Monte 
Carlo (MC) method with importance sampling. In the calculation, 10 6 MC samples are 
taken. 

First, we calculate for six sets of t>2,t>4 and v 6 . The result for the density is displayed in 
Fig. |I|, which shows clearly that the density is not universal. Next, we show the smoothed 
correlations in Figs. |2] and |3|. Again, they are clearly not universal. Now let us multiply 
the smoothed correlation p^ mooth (/i, v) by a 2 and express the result in terms of the scaling 
variables x = \ija and y = vja. The smoothed connected correlation is newly defined by 
the scaling variables as 

pr° oth (x,y) = a 2 p s r oth (P,v) (8) 

Quite dramatically, all the results for p c smooth (x, y) fall on the same universal curve, as shown 
in Figs. £| and ^. 

Brezin and Zee p| have also considered more general probability distribution with even 
polynomials V((f>) and W(<f>) as 

P{(j)) oc exp(-iV Tr V((f>) - (Tr W (0)) 2 ) (9) 

They showed the ensemble with the distribution @ can be mapped to that of the distri- 
bution (|1]) in the large iV limit. Thus the smoothed connected correlation should also obey 
the universal behavior (|5]). 

In this case, we also check the universal behavior numerically. We limit the form of V(<p) 
and W{4>) to 
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V{<t>) = v 2 <p 2 + v^ + v^ 

W{<P) = w 2 <p 2 + + w & <f (10) 

We choose iV to be 100 and take 10 6 MC samples. We calculate for three sets of these 
parameters. The results are shown in Figs. [] and |7|. 

We take smoothing intervals 8 fx and 8v to be around 0.25, which satisfy the smoothing 
conditions mentioned above. Then the dominant factor of errors is due to the MC sampling. 
The errors are large near the endpoint of the spectra. The errors grow as \\x — v\ goes large, 
because the number of the eigenvalues between \x and v increases. There are two regions 
where the theoretical result (H) is not applicable. One is near the endpoint of the spectra 
of the eigenvalues and the other is where two eigenvalues fi and v are close. Except in 
the two regions, the numerical calculations for the distribution ([!]) agree very well with the 
theoretical correlation (|5|) as shown in Figs. f| and [5| The results for the distribution ([|) are 
not so good as those for the distribution ([!]) but consistent with (|) (Figs. |6] and 0). 

The numerical calculations show that the smoothed connected correlation is universal 
when the eigenvalues are rescaled by the endpoint of the spectra of the eigenvalues, although 
the density of the eigenvalues and the bare correlation of the eigenvalues are not universal. 
Thus we conclude that the numerical results agree well with the universal behavior of the 
correlation function @ by Brezin and Zee. Thus the ansatz used in their derivation is well 
established within the numerical errors. 

Morita et. al. @ showed that the density of the eigenvalues are separated into two parts 
in some parameter region. Moreover within one part, the connected correlations of the 
eigenvalues was shown to obey the universal behavior. The numerical calculations in this 
parameter region agree with the theoretical prediction to the extent that they agreed in the 
previous non-separated cases. 

There are also discussions that the correlations are universal for the Wigner ensemble. 
H The Wigner ensemble || is not unitary invariant and each element of the matrix is 
independent and takes 1 or -1. Our numerical calculations for the Wigner ensemble do not 
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agree well with the correlation function (|5|). 

It is a pleasure to thank Y. Avishai for useful discussions. 
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FIGURES 

FIG. 1. The density of the eigenvalues for the probability distribution (1). The param- 
eters (v 2 , «4, ve) are :(a) (0.46, 0.25, 0.21), (b) (0.87, 0.36, 0.68), (c) (-0.70, -0.89, 0.91), 
(d) (5.77, -5.37, 5.62), (e) (3.32, -3.72, 5.20) and (f) (0.0271, -4.63, 2.55). 

FIG. 2. The smoothed connected correlation (fi — u) 2 p s " looth (p, ^)| M =o- The parameters 
(t>2, vq) are the same as in Fig. 1. 

FIG. 3. The smoothed connected correlation (fi — i/) 2 /9^ mooth (^, ^)| M =o.5- The parameters 
(t>2, i»4, vq) are the same as in Fig. 1. 

FIG. 4. The smoothed connected correlation (x — y) 2 p c smooth (x, y)\ x =o, where x and y are the 
scaling variables. The parameters (t>2, V4, vq) are the same as in Fig. 1. Solid line is the theoretical 
result (5). 

FIG. 5. The smoothed connected correlation (x — y) 2 /5 c smooth (x, y)| x =o.5j where x and y are the 
scaling variables. The parameters (t>2, V4, vq) are the same as in Fig. 1. Solid line is the theoretical 
result (5). 

FIG. 6. The smoothed connected correlation (x — y) 2 /5 c smooth (x, y)\ x =o for the probability dis- 
tribution (9), where x and y are the scaling variables. The parameters (t>2, V4, vq, W2, w&, wq) 
are :(a) (1.23, 7.90, 0.290, -0.385, 2.57, 2.80), (b) (3.19, -5.35, 5.66, -6.28, -9.55, 6.77) and 
(c) (7.45, 5.67, 8.94, 3.14, -0.0231, 8.75). Solid line is the theoretical result (5). 

FIG. 7. The smoothed connected correlation (x — y) 2 p c smooth (x, y)| x =o.5 f° r probability distri- 
bution (9), where x and y are the scaling variables. The parameters (v2, V4, vq, u>2, W4, wq) are 
the same as in Fig. 6. Solid line is the theoretical result (5). 
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